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ELECTRICAL NERVE STIMULATION 
BASED ON CHANNEL SPECIFIC SAMPLING SEQUENCES 

Priority 

This application claims priority from United States provisional patent 
application serial number 60/150,773 filed August 26, 1999, entitled Concept 
for Electrical Stimulation of the Acoustic Nerve Based on Channel Specific Sampling 
Sequences (CSSS), the disclosure of which is incorporated herein by reference. 

Field of the Invention 

The present invention relates to electrical nerve stimulation, and more 
particularly, electrostimulation of the nerve based on channel specific 
sampling sequences. 

Background 

Cochlear implants (inner ear prostheses) are a possibility to help 
profoundly deaf or severely hearing impaired persons. Unlike conventional 
hearing aids, which just apply an amplified and modified sound signal, a 
cochlear implant is based on direct electrical stimulation of the acoustic nerve. 
The intention of a cochlear implant is to stimulate nervous structures in the 
inner ear electrically in such a way that hearing impressions most similar to 
normal hearing are obtained. 

A cochlear prosthesis essentially consists of two parts, the speech 
processor and the implanted stimulator. The speech processor contains the 
power supply (batteries) of the overall system and is used to perform signal 
processing of the acoustic signal to extract the stimulation parameters. The 
stimulator generates the stimulation patterns and conducts them to the 
nervous tissue by means of an electrode array which usually is positioned in 
the scala tympani in the inner ear. The connection between speech processor 
and stimulator is established either by means of a radio frequency link 
(transcutaneous) or by means of a plug in the skin (percutaneous). 

At present, the most successful stimulation strategy is the so called 
"continuous-inter leaved-sampling strategy" (CIS), as described by Wilson B. 
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S., Finley C. C, Lawson D. T., Wolford R. D., Eddington D. K., Rabinowitz W. 
M., "Better speech recognition with cochlear implants," Nature, vol. 352, 236 - 
238 (July 1991) [hereinafter Wilson et aL, 1991], which is incorporated herein 
by reference. Signal processing for CIS in the speech processor involves the 
5 following steps: 

(1) splitting up of the audio frequency range into spectral bands by means 
of a filter bank, 

(2) envelope detection of each filter output signal, 

(3) instantaneous nonlinear compression of the envelope signal (map law). 
10 According to the tonotopic organization of the cochlea, each 

stimulation electrode in the scala tympani is associated with a band pass filter 
of the external filter bank. For stimulation, symmetrical biphasic current 
pulses are applied. The amplitudes of the stimulation pulses are directly 
obtained from the compressed envelope signals (step (3) of above). These 

15 signals are sampled sequentially, and the stimulation pulses are applied in a 
strictly non-overlapping sequence. Thus, as a typical CIS-feature, only one 
stimulation channel is active at one time. The overall stimulation rate is 
comparatively high. For example, assuming an overall stimulation rate of 
18kpps, and using an 12 channel filter bank, the stimulation rate per channel 

20 is 1.5kpps. Such a stimulation rate per channel usually is sufficient for 
adequate temporal representation of the envelope signal. 

The maximum overall stimulation rate is limited by the minimum 
phase duration per pulse. The phase duration cannot be chosen arbitrarily 
short, because the shorter the pulses, the higher the current amplitudes have 

25 to be to elicit action potentials in neurons, and current amplitudes are limited 
for various practical reasons. For an overall stimulation rate of 18kpps, the 
phase duration is 27|us, which is at the lower limit. 

Each output of the CIS band pass filters can roughly be regarded as a 
sinusoid at the center frequency of the band pass filter, which is modulated by 

30 the envelope signal. This is due to the quality factor Q ~ 3 of the filters. In 

case of a voiced speech segment, this envelope is approximately periodic, and 
the repetition rate is equal to the pitch frequency. 
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In the current CIS-strategy, the envelope signals only are used for 
further processing, i.e., they contain the entire stimulation information. For 
each channel, the envelope is represented as a sequence of biphasic pulses at 
constant repetition rate. As a characteristic feature of CIS, this repetition rate 
5 (typically l.Skpps) is equal for all channels, and there is no relation to the 

center frequencies of the individual channels. It is intended that the repetition 
rate is not a temporal cue for the patient, i.e., it should be sufficiently high, so 
that the patient does not percept tones with a frequency equal to the 
repetition rate. The repetition rate is usually chosen greater than at twice the 

10 bandwidth of the envelope signals (Nyquist theorem). 

Summary of the Invention 
In accordance with one aspect of the invention, electrodes in a 
multichannel electrode array are activated using channel specific sampling 
sequences. A channel specific sampling sequence for each electrode is 

15 defined, having a particular duration, amplitude, and number of pulses. A 
weighting factor is applied to the channel specific sampling sequence, 
creating a weighted channel specific sampling sequence. Each electrode in 
the multichannel electrode array is then simultaneously activated using sign- 
correlated pulses, the sign-correlated pulses based on parameters of spatial 

20 channel interaction, non-linear compression, and each electrode's weighted 
channel specific sampling sequence. 

In accordance with other related emodiments, the electrodes stimulate 
the acoustic nerve. The multichannel electrode array can be used in a 
monopolar electrode configuration having a remote ground. The pulse 

25 amplitudes can be derived by sampling a signal waveform, for example, one 
half the period of a sinusoid between 0 and n, or one quarter of a sinusoid 
between 0 and k/2 so that the amplitude distribution is monotonically 
increasing. Symmetrical biphasic current pulses can be used to sample the 
signal waveform. The channel specific sampling sequence pulse rate may be 

30 between 5-10 kpps. The parameters of spatial channel interaction can be 

based on a single electrode model having exponential decays of the potentials 
at both sides of the electrode, the sign-correlated pulses having amplitudes 
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which are calculated using properties of a tri-diagonal matrix. The 
multichannel electrode array can be in a cochlear implant, whereby the 
weighting factor is transmitted to the cochlear implant. Start and stop bits, 
and addresses associated with an electrode can also be transmitted to the 
5 cochlear implant. 

In accordance with another embodiment of the invention, electrodes in 
a multichannel electrode array are activated using channel specific sampling 
sequences by applying an acoustic signal to a bank of filters, each filter in the 
bank of filters associated with a channel having an electrode. A weighting 

10 factor is derived for each channel based on the output of each channel's filter. 
The weighting factor is then applied to a channel specific sampling sequence 
having a particular duration, amplitude and number of pulses, creating a 
weighted channel specific sampling sequence. Each channel's electrode is 
simultaneously activated using sign-correlated pulses, the sign-correlated 

15 pulses based on the weighted channel specific sampling sequence, non- 
linear compression, and parameters of spatial channel interaction. 

In accordance with other related embodiments, the electrodes can 
stimulate the acoustic nerve. The weighting factor can be derived by 
rectifying the output of each filter, and then determining the maximum 

20 amplitude of each half-wave in the rectified signal. The multichannel 
electrode array can used a monopolar electrode configuration having a 
remote ground. The pulse amplitudes of the channel specific sampling 
sequence can be derived by sampling a signal waveform, such as one half the 
period of a sinusoid between 0 and n, or one quarter of a sinusoid so that the 

25 amplitude distribution is monotonically increasing. Symmetrical biphasic 
current pulses can be used to sample the waveform. Each channel filter can 
be a bandpass filter. The duration and number of pulses in the channel 
specific sampling sequence can then be derived from the center frequency of 
the channel's bandpass filter. For example, the duration of the channel 

30 specific sampling sequence can be one half of the period of the bandpass 

filter's center frequency. The parameters of spatial channel interaction can be 
based on a single electrode model having exponential decays of the potentials 
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at both sides of the electrode, the sign-correlated pulses having amplitudes 
determined by using properties of a tri-diagonal matrix. The multichannel 
electrode array can be in a cochlear implant, whereby the weighting factor is 
transmitted to the cochlear implant. Start and stop bits, and addresses 
5 associated with an electrode can also be transmitted to the cochlear implant. 

In accordance with another embodiment of the invention, electrodes 
are simultaneously activated in a multichannel electrode array using channel 
specific sampling sequences. Sign-correlated pulses are used. The 
amplitudes of the sign-correlated pulses are calculated by taking into account 
10 parameters of spatial channel interaction. In calculating the amplitudes of the 
sign-correlated pulses a single electrode model having exponential decays of 
the potentials at both sides of the electrode can be used. The amplitudes of 
the sign-correlated pulses can be calculated using properties of a tri-diagonal 
matrix. 

15 In accordance with another embodiment of the invention, channel 

specific sampling sequence having a pulse descriptive characterization are 
defined. The channel specific sampling sequence is used to activate electrodes 
in a multichannel electrode array, each filter in a bank of filters associated 
with a channel having an electrode. Pulse amplitudes of the channel sampling 

20 sequence are derived by sampling a signal waveform. The duration and 

number of pulses of the channel specific sampling sequence are derived from 
a frequency associated with the channel's filter. 

In accordance with other related embodiments, the sampling is of a 
half period of a sinusoid between 0 and tl The sampling can also be of a 

25 quarter period of a sinusoid between 0 and re/ 2, so that pulse amplitude 

distribution monotonically increases. The sampling can use biphasic current 
pulses. Each filter can be a bandpass filter. The duration and number of 
pulses in the channel specific sampling sequence can be derived from the 
center frequency of the charmers bandpass filter. The duration of the channel 

30 specific sampling sequence can be one half of the period of the bandpass 
filter's center frequency. 
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In another embodiment of the invention, a weighting factor for a 
channel specific sampling sequence is derived, the channel specific sampling 
sequence being used to activate electrodes in a multichannel electrode array, 
each filter in a bank of filters associated with a channel having an electrode. 
5 The output of each filter is rectified, creating a half-wave rectified signal. The 
maximum amplitude of each half-wave in the half-wave rectified signal is 
then determined. 

Brief Description of the Drawings 
The foregoing features of the invention will be more readily 
10 understood by reference to the following detailed description, taken with 
reference to the accompanying drawings, in which: 

Figure 1 shows channel specific sampling sequences (CSSS) for two 6-channel 
systems utilizing biphasic pulses at lOkpp/s and phase duration of 25|LLs 
a. derived from a sinusoid within [0 n] 
15 b. derived from a sinusoid within [0 7t/2], amplitudes 

monotonically increasing 

Figure 2 shows stimulation with channel specific sampling sequences (CSSS) 

a. Band pass filter output signal (653Hz - 876Hz) 

b. Half wave rectified band pass filter output 
20 c. Associated CSSS stimulation sequence 

Figure 3 shows stimulation with channel specific sampling sequences (CSSS) 

a. Bandpass filter output signal (3457Hz -5500Hz) 

b. Half wave rectified band pass filter output 

c. Associated CSSS stimulation sequence 

25 Figure 4 shows a comparison of stimulation patterns between CSSS and CIS 

a. Band pass filter output signal (653Hz - 876Hz) 

b. CSSS stimulation sequence 

c. CIS stimulation sequence (envelope sampling with single pulses 
at 1.5 kpps) 

30 Figure 5 shows a comparison of stimulation patterns between CSSS and CIS 

a. Band pass filter output signal (3457Hz - 5500Hz) 

b. CSSS stimulation sequence 



6 



WO 01/13991 PCT/IB00/01338 

c. CIS stimulation sequence (envelope sampling with single pulses 
at 1.5 kpps) 

Figure 6 shows estimated potential distributions in the scala tympani (A, = 
3.6mm, d = 2.8mm) 
5 a. Responses to single channel activation 

b. Effective Potential Distributions (asterisks for CIS, circles for 

CSSS) 

Figure 7 shows a comparison of overall stimulation patterns between CSSS 
and CIS (electrode distance: d = 2.8mm, space constant: X =3.6 mm) 

10 a. 6-channel CSSS 

b. 6-channel CIS 

Detailed Description of the Invention 
A cochlear implant with stimulation patterns containing enhanced 
temporal information, especially in the low frequency range up to 1kHz, is 

15 described. It is known from literature that the neurons are able to track 

analogue electrical sinusoidals up to about 1kHz. This ability is not exploited 
in the present CIS strategy, since the sampling rate is too low to represent 
high frequency envelope waveforms. 

The stimulation strategy utilized is based on channel specific sampling 

20 sequences (CSSS). The basic idea is to a apply a stimulation pattern, where a 
particular relationship to the center frequencies of the filter channels is 
preserved, i.e., the center frequencies are represented in the temporal 
waveforms of the stimulation patterns, and are not fully removed, as in CIS. 

Each stimulation channel is associated with a particular CSSS, which is 

25 a sequence of ultra-high-rate biphasic pulses (typically 5-10kpps). Each CSSS 
has a distinct length (number of pulses) and distinct amplitude distribution. 
The length of a CSSS is derived from the center frequency of the associated 
band pass filter. A CSSS associated with a lower filter channel is longer than a 
CSSS associated with a higher filter channel. Typically, it is one half of the 

30 period of the center frequency. The amplitude distribution can be adjusted to 
patient specific requirements. For convenience, the amplitude of the 
maximum biphasic pulse within a CSSS is normalized to one. For illustration, 
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two examples for a 6-channel system are shown. In Fig. 1(a), the CSSS's are 
derived by sampling one half of a period of a sinusoid, whose frequency is 
equal to the center frequency of the band pass filter (center frequencies at 
440Hz, 696Hz, 1103Hz, 1745Hz, 2762Hz, and 4372Hz). Sampling is achieved 
5 by means of biphasic pulses at a rate of lOkpps and a phase duration of 25]tis. 
For channels #5 and #6, one half of a period of the center frequencies is too 
short to give space for more than one stimulation pulse, i.e., the "sequences" 
consist of only one pulse, respectively. In Fig. 1(b), the sequences are derived 
by sampling one quarter of a sinusoid with a frequency, which is half the 

10 center frequency of the band pass filters. These CSSS's have about the same 
durations as the CSSS's in Fig. 1(a), respectively, but the amplitude 
distribution is monotonically increasing. Such monotonic distributions might 
be advantageous, because each pulse of the sequence can theoretically 
stimulate neurons at sites which cannot be reached by its predecessors. This 

15 is a pure "geometric" effect, and could possibly result in a broader temporal 
distribution of the firing pattern of the neurons. 

An example of a stimulation pattern based on CSSS is depicted in Fig.2 
for a voiced speech segment. For reasons of clarity, the influence of spatial 
channel interaction is neglected here. In addition, and in the following text, 

20 the instantaneous non-linear compression is omitted for convenience, 

however it is realized that such conversion is required for actual stimulation 
patterns. Fig. 2(a) shows the output of a band pass filter (cut off frequencies 
at 553Hz and 876Hz). Fig. 2(b) illustrates the half-wave rectified version of the 
signal. In Fig. 2(c), each half-wave-pulse is replaced by a CSSS, where the 

25 amplitude of the maximum pulse within each CSSS is equal to the maximum 
of the associated half-wave-pulse. Thus, Fig. 3 represents a sequence of 
weighted and time-shifted CSSS's. The CSSS used for this example is equal to 
the CSSS in Fig. 1(a) for channel CH2, and for convenience, each biphasic 
pulse is represented as a single vertical line. 

30 An example of a stimulation pattern based on CSSS for a higher 

frequency channel is shown in Fig. 3 (the input speech segment is the same as 
for Fig. 2, spatial channel interaction is neglected again). The band pass filter 
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here selects a range between 3475Hz and 5500Hz. With a center frequency of 
4273Hz, the period is 229|^s, and sampling one half of this period gives space 
for only one pulse (cf. CSSS as shown in Fig. 1 for channel CHS). Here, the 
envelope sampling is reduced to a sampling with single pulses at a rate equal 
5 to about the center frequency of 4273Hz. 

In Fig. 4 stimulation sequences of the new approach are directly 
compared to the corresponding CIS-sequences at l.Skpps. The CSSS-based 
sequence in Fig. 4(b) clearly represents the temporal fine structure plus the 
envelope information of the band pass output shown in Fig. 4(a), whereas the 

10 CIS-pattern in Fig. 4(c) is obtained by sampling the envelope, and thus any 
temporal fine structure is removed. At a stimulation channel at higher 
frequencies, Figs. 5(b) and (c) are derived by envelope sampling with single 
pulses. However, in this frequency range, neurons are only able to track the 
envelope signals, but cannot follow the stimulation frequency itself. The 

15 difference between traces 2 and 3 is the sampling rate, which is considerably 
lower for CIS. 

For the practical implementation of the new stimulation approach as 
described above it is necessary to utilize simultaneous stimulation techniques. 
Interleaved sampling as employed in CIS is impractical here, since this would 
20 require stimulation pulses with a phase duration of only few microseconds. 

In the following, the most important mechanisms of channel 
interaction involved with cochlear implants are summarized. 

Channel interaction 

In principle, channel interaction in connection with pulsatile 
25 stimulation strategies occurs as a spatial and as a temporal effect. Temporal 
interaction could be further separated into "physical" and "physiological" 
interaction. 

a. Spatial channel interaction 

Spatial channel interaction means that there is considerable geometric 
30 overlapping of electrical fields at the location of the excitable nervous tissue, if 
different stimulation electrodes (positioned in the scala tympani) are 
activated. Thus, neglecting temporal channel interaction, the same neurons 
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can be activated if different electrodes are stimulated. Stimulation of a 
particular electrode against a remote ground electrode (monopolar 
stimulation) causes an electrical potential within the scala tympani which can 
roughly be described by two decaying exponentials at both sides of the 
5 electrode, and the space constant (in humans) is typically X = 3.6mm, as 

described by Wilson B. S., Finley C. C, Zerbi M., and Lawson D. T., "Speech 
processors for auditory prostheses/' Seventh Quarterly Progress Report, 
Feb.l st through April 30 th , 1994, NIH Contract N01-DC-2-2401 [hereinafter 
Wilson et al., 1994], which is incorporated herein by reference. This type of 

10 channel interaction is first of all due to the conductive fluids and tissues 

surrounding the stimulation electrode array. A similar space constant is also 
obtained by simulation, if a simple model of a cochlea composed of 
exclusively ohmic resistors is assumed, as described by Krai A., Hartmann R., 
Mortazavi D., and Klinke R., "Spatial resolution of cochlear implants: the 

15 electrical field and excitation of auditory afferents," Hearing Research 121, 
pp. 11-28, (1998), which is incorporated herein by reference. This model 
allows a rough quantitative computation of the electrical potentials within the 
scala tympani, as well as at the position of excitable neurons, 
b. Physical temporal channel interaction 

20 Physical temporal channel interaction means that the electrical 

properties of a stimulation pulse in the nervous tissue are biased by its 
predecessor, e.g., due to residual charge stored in the tissue and in the 
membrane capacitances of the neurons. Physical temporal interaction is 
suppressed to a great extent by using symmetrical, biphasic stimulation 

25 pulses. Most of the charge delivered to the excitable tissue during the first 
phase of a stimulation pulse is removed during the second. However, since 
the tissue shows some capacitative behavior, some residual charge remains 
after the end of the stimulation pulse and possibly may bias the subsequent 
stimulation pulse. Theoretically, triphasic pulses (with zero net charge) 

30 would help to further reduce physical temporal channel interaction. 
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c. Physiological temporal channel interaction 

Physiological interaction means effects associated with the refractory 

properties of the neurons. Following Wilson et al, 1994, a recovery function 

r(t) can be defined as 
5 r(t) = 0, for t < a, and 

r(t) = l-exp(--^H fort>t a/ (1) 
x 

with an absolute refractory period t a ~ 700(is, and a time constant 
x ~ 1250(is for the relative refractory period. For example, if two supra- 
threshold stimulation pulses are applied, and the second pulse falls into the 

10 absolute refractory period after the first, no additional action potential can be 
elicited. If the second pulse occurs during the relative refractory period, an 
enhanced amplitude is necessary to generate an action potential. 

The influence of physiological temporal interaction on speech 
understanding is currently investigated at various research centers 

15 worldwide. At the moment, it seems that the similarity between neural 

excitation patterns due to electrical stimulation and natural excitation patterns 
can be enhanced, if very high stimulation rates are employed (> 3kpps per 
channel, as described by Matsuoka A. J., ''Compound action potentials 
evoked by electrical pulse trains: effects of stimulus parameters on response 

20 patterns/ 7 thesis at University of Iowa, Quly 1998), which is incorporated 
herein by reference. High rates may mimic membrane noise (spontaneous 
activity) and thereby keep different neurons in different refractory states. If 
this is the case, it can be expected that the ensemble spiking patterns can 
reflect the envelope of amplitude modulated electrical pulse sequences up to 

25 considerably higher frequencies, and thus more temporal information can be 
provided to the brain. 

Consideration of spatial channel interaction 

In CIS strategy, the influence of spatial channel interaction is reduced 
by employing pulses which are not overlapping in time (interleaved 
30 sampling). The conductivity in the scala tympani here leads to a considerable 
spread and a de-focusing of the electrical field at the site of the excitable 
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tissue. However, an additional effect occurs, if simultaneous stimulation of 
two or more electrodes against a remote ground electrode is considered. Here 
the conductivity represents a shunt conductance between active electrodes, 
which in general results in a temporal mixture of constructive and destructive 
5 superposition of electrical fields at the position of the neurons. For example, 
if two simultaneous stimulation channels produce currents with equal 
amplitudes, but different signs, most of the current will flow through the 
shunt conductance and will not reach the intended neurons. This additional 
effect can be removed, if "sign-correlated" pulses are employed. 
10 Sign-correlation here means that the signs of the phases of simultaneous 

stimulation pulses are equal. This ensures that the sum of the magnitudes of 
the single stimulation currents is forced to flow into the reference electrode. 
Thus, at the site of the excitable neurons, only constructive superposition of 
currents is possible. 

15 The injection of a current by means of a single active electrode into the 

scala tympani causes a particular voltage in the tissue just close to the 
electrode (measured against the remote reference electrode), and an 
exponential decay at both sides of the electrode. The space constant typically 
is A, = 3.6mm, as described by Wilson et al, 1994. Assuming a linear and pure 
20 ohmic system, the injection of currents in more than one electrode causes a 
superposition of the potential distributions due to the single currents. 

The idea here is to modify stimulation currents such that at least the 
potentials at the position of the electrodes are equal as in the case of single 
channel stimulation. Assuming N channels, the single channel (non- 
25 simultaneous) current amplitudes x n (n = 1-N) and the amplitudes y n (n = 1-N) 
for simultaneous channels are related via the following set of linear equations: 
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where H 1 is the inverse matrix of H. Fortunately, matrix H 1 in general 
is a tri-diagonal matrix with non-zero elements only in the main-, the upper 
and lower neighboring diagonals (see Appendix). 

An example is shown in Fig.6 for six electrodes (N = 6). The x-axis is 
normalized to a distance d = 2.8mm between the electrodes, i.e., the electrodes 
are at positions 1 to 6. A space constant X = 3.6mm is assumed. The y-axis is 
normalized to the maximum potential of electrode #4 at position 4. Fig.6(a) 
depicts the single voltage distributions in the scala tympani as responses to 
single electrode currents at different amplitudes. 

For CIS, the electrodes are activated sequentially, and thus each of the 
single potential distribution applies for the duration of a pulse phase. 
Assuming a pulse repetition rate of 1.5kppulses/s for each channel, the 
overall time necessary to present all six distributions is 666|is, which is just 
about the duration of the absolute refractory period (t a ~ 700|~is). This allows 
the following rough approximation: for CIS, due to physiological channel 
interaction, the "effective" stimulation pattern is the contour of the single 
potential distributions, as shown in Fig.6(b), (asterisks). 

For Fig.6(b) (circles), the amplitudes y n (n = 1-6) for simultaneous 
stimulation are computed by means of (4). As demanded, the potentials 
coincide at the electrode positions. Obviously, the peaks obtained by taking 
the contour of the non-simultaneous potential distributions CIS are more 
pronounced than with CSSS. Unfortunately, not all amplitude distributions 
x n > 0 yield solutions y n with positive elements for all n. This is in 
contradiction to the principle of "sign-correlation", and requires to compute a 
modified vector y' n , which contains only non-negative elements (see 
Appendix). 

For CIS, the electrodes are activated sequentially, and thus each of the 
single potential distribution applies for the duration of a pulse phase. 
Assuming a pulse repetition rate of 1.5kppulses/s for each channel, the 
overall time necessary to present all six distributions is 666|lls, which is just 
about the duration of the absolute refractory period (a ~ 700|Lis). This allows 
the following rough approximation: for CIS, due to physiological channel 
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interaction, the "effective" stimulation pattern is the contour of the single 
potential distributions, as shown in Fig. 6(b), (asterisks). 

If amplitudes I iCSSS (i = 1-6) are computed by means of (4) for a 
distribution of I (i = 1-6) as in Fig.6, one of the current amplitude, I 5 CSSS , is 
5 negative. This is in contradiction to the concept of sign-correlated pulses, 
which requires the same sign for all simultaneously activated channels. In 
this case, e.g., the 5 th row and 5 th column of Matrix H are deleted, resulting in a 
modified 5x5-Matrix H mod . A distribution I iCSSSmod with only 5 simultaneous 
currents is obtained by (1) computing the inverse matrix H _1 mod , and (2) 

10 multiplying H' 1 ^ with vector I (i = 1,2,3,4, and 6). The resulting potential 
distribution is shown in Fig. 6(b), (circles). As demanded, the potentials 
coincide at the positions 1, 2, 3, 4, and 6. Obviously, the peaks obtained by 
taking the contour of the non-simultaneous potential distributions CIS are 
more pronounced than with CSSS. 

15 For the CSSS system, Fig. 7(a), the envelope sampling sequences for 

each channel are chosen as shown in Fig. 1(a). To obtain the actual stimulation 
signals for each channel, the spatial channel interaction is taken into account. 
As expected, the stimulation pattern reflects the temporal fine structure. In 
particular, the center frequency of channel #2 is represented in the temporal 

20 waveform. A so-called "hole-effect" can be observed: if electrode #2 is not 
active, i.e., if the output of filter channel #2 is negative, then other spectral 
maxima are not masked (due to spatial channel interaction) and appear in the 
waveform. 

The CIS system, Fig. 7(b) is based on an overall sampling rate of 
25 lOkpps, resulting in a rate of 1667pps per channel for the 6-channel system. 
Envelope detection for each channel is achieved with a full-wave rectifier and 
a low pass filter with a cut off frequency of 400Hz (Butterworth filter of 2 nd 
order), respectively. Obviously, the envelope signal is sampled and presented, 
but the temporal fine structure is lost. 
30 Implementation of a cochlear implant system based on CSSS 

Although based on highly synchronous stimulation, the CSSS 
approach is well suited for an implementation in a practical cochlear implant 
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system. The information transfer rate from the speech processor to the 
implant can be kept comparatively low. An example of a data word for a 
12-channel CSSS system is shown in Tab.l. 

One data word consists of 16 bits, including START and STOP bits. The 

5 two special bits SPEC1 and SPECO represent a low rate information channel 
and are used for the initialization of the implant. The implant is permanently 
supplied with data defining the amplitude distributions and pules repetition 
rate of the normalized CSSS's, as well as data defining the reference current 
levels for each channel. These data are stored in particular implant memories. 

10 Besides, safety bits (e.g., for cyclic redundancy check (CRC)) are transmitted. 
Note that for proper operation of the implant, the information defining the 
normalized 

CSSS's and the reference current levels theoretically have to be transmitted 
only once. 

15 The four address bits ADD3-ADD0 define the channel address, and 

bits W7-W0 the weighting factor of the associated CSSS. The repetition rate of 
the 



Tab.l: Data word 



Bit# 


Definition 


1 


START 


2 


ADD3 


3 


ADD2 


4 


ADD1 


5 


ADDO 


6 


W7 


7 


W6 


8 


W5 


9 


W4 


10 


W3 


11 


W2 


12 


Wl 



16 



WO 01/13991 



PCT/IB00/01338 



13 


WO 


14 


SPEC1 


15 


SPECO 


16 


STOP 



CSSS's is comparatively low, especially at the low frequency channels. It is 
not necessary to transmit the amplitudes of the individual ultra-high rate 
pulses, since the amplitude distributions are already stored in the implant. 
5 Assuming an input analysis range between 350Hz - 5500Hz for a 12- 

channel system, and a logarithmic spacing of the band pass filter ranges, 
results in center frequencies 393Hz, 494Hz, 622Hz, 782Hz, 983Hz, 1237Hz, 
1556Hz, 1958Hz, 2463Hz, 3098Hz, 3898Hz, and 49036Hz. Thus, the average 
CSSS-repetition rate is equal to the sum of the center frequencies, i.e., 

10 R csss = 22386Hz. This is equal to the average data word repetition rate R dataword . 
The resulting average overall bit rate is R bit = 16R datavvord ~ 358kbit/s. Thus, a bit 
rate of 600kbit/s for a practical cochlear implant is sufficient for complete 
information transfer. However, this is a moderate rate as compared to the 
case, if each stimulation pulse has to be defined independently. Here, 

15 assuming a frame-rate of lOkpps of simultaneous stimulation pulses and a 
data word of 16bit per pulse, an overall bit rate of 1920kbit/s results. Such a 
bit rate is almost impossible to realize with a inductive link system at 
reasonable power consumption. 

Within the implant, the correction of the amplitudes due to spatial 

20 channel interaction has to be performed for each simultaneous stimulation 
frame. 

Summary 

In summary, the CSSS stimulation approach may be summarized as 
follows. 

25 (1) For stimulation, a multichannel electrode array within the scala 

tympani and a remote ground electrode is used (monopolar electrode 
configuration). The basic stimulation waveform is a symmetrical, biphasic 
pulse. 
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(2) Stimulation involves simultaneous activation of electrodes in the 
scala tympani employing sign-correlated pulses. Sign-correlated means that if 
two or more pulses occur simultaneously at different electrodes, positive and 
negative phases are absolute synchronous in time. 
5 (3) The amplitudes of the sign-correlated pulses are estimated by 

taking into account parameters of spatial channel interaction. Assuming that 
a single electrode causes exponential decays of the potentials at both sides of 
the electrode allows a computationally efficient calculation of the pulse 
amplitudes, since a tri-diagonal matrix is involved. 
10 (4) Processing of the acoustic signal involves a filter bank for 

splitting up the audio frequency range (similar to CIS). According to the 
tonotopic organization of the scala tympani, each band pass filter is associated 
with a stimulation electrode. 

(5) Each stimulation channel is associated with a normalized, 
15 channel specific sampling sequence (CSSS) of ultra-high-rate pulses. 

Typically, rates between 5-10kpps are employed. For each channel, the CSSS 
has different length and different amplitude distribution. The maximum 
amplitude of a normalized CSSS is one. 

(6) The length of a CSSS is derived from the center frequency of the 
20 associated band pass filter. Typically, it is one half of the period of the center 

frequency. For example, a band pass center frequency of 500Hz results in a 
CSSS-length of 1ms comprising 10 pulses, (assuming a ultra-high-rate of 
lOkpps). 

(7) The amplitude distribution of a CSSS is chosen for optimum 

25 performance with respect to mimicking membrane noise. As many neurons as 
possible shall be kept in different refractory states. 

Although various exemplary embodiment of the invention have 
been disclosed, it should be apparent to those skilled in the art that various 
changes and modifications can be made which will achieve some of the 
30 advantages of the invention without departing from the true scope of the 

invention. These and other obvious modifications are intended to be 
covered by the claims that follow. 
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The matrix product (2) can be regarded as a convolution product of an 
infinite sequence h n and a sequence y n , with non-zero elements only at 
5 positions n = 1, 2, ...N, i.e., 

(Al) 

10 where sequence h n is given by 

h n =a n u n +cT n u_ n _ 1 . 

(A2) 

15 Function u n is the unit step, i.e., u n = 1 for n > 0, and u n = 0 for n < 0. 

Sequence h n represents an infinite impulse response (IIR) with exponential 
decays at both sides of the origin ( I a I < 1). The z-transform is given by 

H(z)= + -^—p, 

(1-oz -1 ) (1-cT'z -1 ) 

20 (A3) 

which can be expressed as 



(« + -) 
H(z) = a 



(z +1 -(a+^ + z -1 ) 



a 

25 (A4) 



Transformation of (Al) into the z-domain yields 
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and thus 

5 

Y(z) = H-'(z)X(z). 

(A6) 

Inserting (A4) yields 

10 

Y(z) = l -r(z +i - (« + -) + z" 1 )X(z) . 

cx 

(A7) 

The inverse z-transform immediately yields 

15 

y n = ^(S n+ , -«x + -)8 n + 5 n _,)*x„, 

(«+-) 
a 

(A8) 

where 6 n is the unit impulse, i.e., 5 n = 1 for n = 0, and 5 n = 0 elsewhere. 
The first term of the convolution product (A8) is a finite impulse response 
20 (FIR). Equation (A8) can be expressed as 

1 / / U 

y n = j~ (x "+' "(a+-)x„ +x n _,), 

a 

(A9) 

25 which is a set of linear equations. To calculate y n at positions n = 1 and 

n = N requires to know amplitudes x 0 and x N+1 . Since sequence y n can have 
non-zero elements only at positions n = 1, 2, ...N, it follows with (Al) 
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x 0 =y,cc+y 2 a 2 +... + y N a N =a(y, +y 2 a' +... + y N a N ') = ax,, 

(A10) 

and similarly 

x n+i =yi aN +y2« N_1 +... + y N a = a(y,a N "' +y 2 a N_2 +... + y N ) = ax N 
(All) 

Inserting x 0 and x N+1 in (A9) for n = 1 and n = N allows to write (A9) as 
matrix equation, and the matrix necessarily has to be identical to the inverse 
matrix of H: 



r y. 






r x, 


\ 


y 2 






x 2 




y 3 




= H 1 






y N - 2 






X N-2 




y N -i 






X N-1 




v yN 


J 




v X N 


/ 



(A12) 

where matrix H' 1 is a tri-diagonal matrix given by 



bo 


— a 


0 .. 


. 0 


0 


0 


— a 


b 


— a .. 


. 0 


0 


0 


0 


-a 


b .. 


0 


0 


0 


0 


0 


0 .. 


. b 


— a 


0 


0 


0 


0 .. 


. —a 


b 


— a 


0 


0 


0 .. 


. 0 


-a 


b 0 



(A13) 
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b„=- 



1 1 



( a) 

a 



a 



b =— — ! (a + — ), and 

a 



a = - 



( a) 

a 



(A14) 



It shall be mentioned that the analysis can simply be expanded to the 
case, if the infinite sequence h n (A2) is of the form 

10 

h n =a n u n +p- n u_ n _„ 

(A15) 



i.e., the exponential decays are different for n > 0 and n < 0 ( I a I < 1, 
15 Ipl <1). 

In the following, it is assumed that for a given vector x n with x n > 0 for 
all n (n = 1, 2, ... N), equation (3) yields a vector y n containing negative 
elements at positions k, i.e., y k < 0. Negative elements mean negative current 
amplitudes, which are in contradiction to the principle of sign-correlation, 
20 and therefore have to be avoided. 

One method of handling such a case is to compute a new vector y' n , 
where the elements at the positions k are set to zero, i.e., y' n=k = 0. This 
restriction requires a modified input vector x' n . In the proposed method, x' n 
differs from vector x n only at positions k and remains unchanged elsewhere, 
25 i.e., x' n;ek = x n9ftk , and x' n=k ^ x n==k . 

It is claimed that conditions 
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yield a definite solution for vector x' n at all positions. 
5 To prove this claim for an arbitrary pattern of k, "zero-sequences" of 

neighboring indices within k of length L are regarded. For convenience, the 
smallest index within each zero-sequence is designated as starting index kO. 
For example, for N = 12, and assuming a pattern k = [1, 2, 5, 6, 7, 10], three 
zero-sequences [1, 2], [5, 6, 7], and [10] with lengths L = 2, 3, and 1 can be 
10 identified, and the starting indices are 1, 5, and 10, respectively. A zero- 
sequence of length L = 1 is also designated as "sequence". 
Two categories have to be distinguished: 

Category (1): a zero-sequence does not contain indices 1 or N, and 
Category (2): a zero-sequence contains either index 1 or N. 
15 For the example of above, zero-sequence [1, 2] belongs to category (2), 

zero-sequences [5, 6, 7], and [10] belong to category (1). 

ad Category (1): here, for a given zero-sequence, neighboring elements 
with positive y n at the lower- and the upper range of a zero-sequence do exist 
at positions n = k0-l and n = kO+L, respectively. For example, for N = 12 and 
20 a zero-sequence [5, 6, 7], kO = 5, and L = 3, and thus the neighboring positions 
are n = 4 and n = 8. 

Setting y' n=k = 0 yields the following set of equations: 

0 = -ax k0 _, + bx\ Q ~ax\ 0+l 

aX k0+L-3 + & X k0+L-2 ~~ aX k0+L-\ 

~ aX k0+L-2 + *o+L-l ~ aX k0+L ' 

25 

(A17) 



0 = 
0 = 
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Elements x k(M and x k0+L/ and coefficients a and b are known, and thus for 
L > 1, (A17) can be written as 



x kO 








' X k0-1 


\ 


X 'k0+1 








0 




X 'k0+2 




= "~Ql 

a 




0 




X kO+L-3 




-l 


0 




X kO+L-2 








0 




^ X 'k0+L-1 


J 






^ X kO+L 


J 



(A18) 



with matrix square Q L 



Ql = 





f b 


-a 


0 .. 


. 0 


0 


o A 




— a 


b 


— a .. 


. 0 


0 


0 


1 


0 


-a 


b .. 


. 0 


0 


0 


















0 


0 


0 .. 


b 


— a 


0 




0 


0 


0 .. 


. —a 


b 


-a 




1° 


0 


0 .. 


0 


-a 


b J 



(A19) 



The number of lines (and rows) of matrix Q L is L (L > 1). Note that 
amplitudes x' k are fully determined by the "neighboring" amplitudes x k0 _j and 
x k0+L . In particular, amplitudes x' k0 and x' k0+L ^ can be calculated with 



_c (L) d (L) 

X kO - X k0-1 H X kO+L / and 

a a 

_ d (D C (L) 
X kO+L-1 ~" X k0-1 H X kO+L ' 

a a 



(A20) 
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where coefficients c (L) and d (L) are the elements at the left- and right 

upper corner of matrix Q" 1 , respectively, i.e., at matrix positions (1,1) and (1,L). 

For each length L, there exists one unique pair of coefficients c (L) and d (L) . For 

a 2 

L = 1, evaluation of (A17) yields c (1) = c (2) = — . With (A20), the amplitudes 

b 

y' k01 and y' k0+L can be determined: 

y'k()-l - " aX kO-2 + k X kO-l _ ax 'k0 ~ 

= -ax k0 _ 2 +(b-c (L) )x k0 _ 1 -d (L) x k0+L ,and 

y kO+L = " aX kO+L-1 + t >X kQ+L ~ aX kCM-L+I = 

= -d (L) x k0 _, +(b-c (L) )x k0+L — ax kQ+L+1 . 

(A21) 

Thus, setting the amplitudes y' k = 0 for a zero-sequence results in a 
modification of the elements in y n only at positions, which are neighboring to 
the zero-sequence. Note that other elements of y n are not concerned. Equation 
(A21) can be implemented by means of the following steps: 

(1) replace coefficients -a, b, and -a in line kO-1 by -a, b+c (L) , 
and - d (L) , 

(2) replace coefficients -a, b, and -a in line kO+L by - d (L) , 
b+c (L) , and -a, 

(3) delete lines and rows with indices k from matrix H' 1 , and 
remove elements with indices k from vector x n . 

ad case (2) : if a zero-sequence contains index 1, the modified 
amplitudes are exponentials up to index L (cf . (A10)), and can be derived 
from amplitude x L+1 : 
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V" L J \ aX L+l 



(A22) 
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Regarding the matrix operation, the coefficient b of line L+l of H" 1 has 
to be replaced by coefficient b 0 . Then all lines and rows with indices k have to 
be removed, and the elements of x n with indices k can be ignored. 

Similarly, if a zero-sequence contains index N, the modified 
amplitudes are exponentials for indices greater than kO+L and can be derived 
from amplitude x k(M : 



k0 



kO+1 



OCX 



kO-1 



a x k0 _j 



J v 



a x k0 _j 



(A23) 



Regarding the matrix operation, the coefficient b of line kO-1 of H" 1 has 
to be replaced by coefficient b 0 . Then all lines and rows with indices k have to 
be removed, and the elements of x n with indices k can be ignored. 

Theoretically, vector y' n can again contain negative elements, but the 
magnitude of the negative elements are comparatively small. A repetition of 
the proposed procedure could remove them, but in many cases, it is sufficient 
to replace the negative elements by zeros and neglect the impact. 

Based on the analysis of above, the following computational efficient 
procedure for the consideration of channel interaction in an N channel system 
can be applied. 
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(1) Compute y n by multiplication of H 1 and x n . 

(2) Select elements y n=k < 0 and set y k ' = 0. 

(3) Modify elements of H' 1 according to (A21), (A22), (A23) 

(4) Remove all lines and rows of H 1 with indices k, and remove all 
elements x k . 

(5) Compute elements y n ', which are neighboring to zero-sequences. 
Example: 

Let the result of the matrix multiplication y n = H J x n (matrix H 1 defined 
by coefficients b 0 , b, and a, for a 12-channel system (N = 12) be a vector 
containing negative elements at positions k = [1, 2, 6, 7, 9, 10, 11]. Then the 
modified vector y' n is 

( 0 > 
0 

y*3 
y' 4 

y's 
o 

y'n - 0 / 

y's 
o 

o 

0 

,y', 2 j 

(A24) 

and the unknown elements are computed by 
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f y ' 3 ' 


f 


y' 4 




y' 5 




y's 




U'l2 j 


V 



-a 



0 
0 
0 



— a 
b 

— a 
0 
0 



0 

— a 

b-c (2) 
-d (2) 
0 



0 
0 

-d (2) 

b _ c (2)_ c (3) 
_ d <3) 



0 
0 
0 

-d (3) 

b 0 -c <3) 



V 



(A25) 



Note that element y' 4 = y 4V because position n = 4 is not neighboring to 
a zero-sequence. Element y'g is neighboring to two zero-sequences. Therefore, 
the corresponding element in the main diagonal is b-c (2) -c (3) , reflecting the 
influence of both zero-sequences. Coefficients c (2) , d (2) and c (3) , d (3) are computed 
by inverting matrices Q 2/ and Q 3/ which themselves only depend on 
coefficients a and b. 
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1. A method of activating electrodes in a multichannel electrode array 
using channel specific sampling sequences, the method comprising: 

a. defining for each electrode a channel specific sampling sequence 
having a selected duration, amplitude, and number of pulses; 

b. applying a weighting factor to each channel specific sampling 
sequence, creating a weighted channel specific sampling sequence for each 
electrode; and 

c. simultaneously activating each electrode using sign-correlated 
pulses, the sign-correlated pulses being based on: 

i. parameters of spatial channel interaction; 

ii. each electrode's weighted channel specific sampling 
sequence; and 

iii. non-linear compression. 

2. A method according to claim 1, wherein the electrodes stimulate the 
acoustic nerve. 

3. A method according to claim 1, wherein the multichannel electrode 
array uses a monopolar electrode configuration having a remote ground. 

4. A method according to claim 1, wherein the channel specific sampling 
sequence pulse amplitudes are derived by sampling a signal waveform. 

5. A method according to claim 4, wherein the sampling is of a half 
period of a sinusoid between 0 and n. 

6. A method according to claim 4, wherein the sampling is of a quarter 
period of a sinusoid between 0 and tc/2 so that pulse amplitude distribution 
monotonically increases. 
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7. A method according to claim 4, wherein the sampling uses 
symmetrical biphasic current pulses. 

5 8. A method according to claim 1, wherein the channel specific sampling 
sequence has a pulse rate between 5-10 kpps. 

9. A method according to claim 1, wherein the parameters of spatial 
channel interaction are based on a single electrode model having exponential 

10 decays of the potentials at both sides of the electrode. 

10. A method according to claim 9, wherein the sign-correlated pulses 
have amplitudes determined by using properties of a tri-diagonal matrix. 

15 11. A method according to claim 1, wherein the multichannel electrode 
array is in a cochlear implant and the weighting factor is transmitted to the 
cochlear implant. 

12. A method according to claim 11, wherein start and stop bits are 
20 transmitted to the cochlear implant, along with the weighting factor. 

13. A method according to claim 11, wherein addresses associated with an 
electrode are transmitted to the cochlear implant, along with the weighting 
factor. 

25 

14. A method of activating electrodes in a multichannel electrode array 
using channel specific sampling sequences, the method comprising: 

a. applying an acoustic representative electrical signal to a bank of 
filters, each filter in the bank of filters associated with a channel having an 

30 electrode; 

b. deriving a weighting factor for each channel from the output of 
each channel filter; 
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c. applying the weighting factor to a channel specific sampling 
sequence having a selected duration, amplitude and number of pulses, 
creating a weighted channel specific sampling sequence; and 

d. simultaneously activating each channel's electrode using sign- 
correlated pulses, the sign-correlated pulses being based on: 

i. parameters of spatial channel interaction; 

ii. each electrode's weighted channel specific sampling 
sequence; and 

iii. non-linear compression. 

15. A method according to claim 14, wherein the electrodes stimulate the 
acoustic nerve. 

16. A method according to claim 14, wherein deriving the weighting factor 
comprises: 

a. rectifying the output of each filter, creating a half-wave rectified 
signal; and 

b. determining a maximum amplitude of each half-wave in the 
half- wave rectified signal. 

17. A method according to claim 14, wherein the multichannel electrode 
array uses a monopolar electrode configuration having a remote ground. 

18. A method according to claim 14, wherein the pulse amplitudes of the 
channel specific sampling sequence are derived by sampling a signal 
waveform. 

19. A method according to claim 18, wherein the sampling is of a half 
period of a sinusoid between 0 and 7C. 
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20. A method according to claim 18, wherein the sampling is of a quarter 
period of a sinusoid between 0 and rc/2 so that pulse amplitude distribution 
monotonically increases. 

5 21. A method according to claim 18, wherein the sampling uses 
symmetrical biphasic current pulses. 

22. A method according to claim 14, wherein each filter is a bandpass 
filter. 

o 

23. A method according to claim 22, wherein the duration and number of 
pulses in the channel specific sampling sequence is derived from the center 
frequency of the channel's bandpass filter. 



15 24. A method according to claim 23, wherein the duration of the channel 
specific sampling sequence is one half of the period of the bandpass filter's 
center frequency. 

25. A method according to claim 14, wherein the parameters of spatial 

20 channel interaction are based on a single electrode model having exponential 
decays of the potentials at both sides of the electrode. 

26. A method according to claim 25, wherein the sign-correlated pulses 
have amplitudes determined by using properties of a tri-diagonal matrix. 

25 

27. A method according to claim 14, wherein the multichannel electrode 
array is in a cochlear implant, and the weighting factor is transmitted to the 
cochlear implant. 



30 28. A method according to claim 27, wherein start and stop bits are 
transmitted to the cochlear implant along with the weighting factor. 
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29. A method according to claim 27, wherein addresses associated with an 
electrode are transmitted to the cochlear implant along with the weighting 
factor. 

5 30. A method of simultaneously activating electrodes in a multichannel 
electrode array using channel specific sampling sequences, the method 
comprising: 

a. using sign-correlated pulses; and 

b. calculating the amplitudes of the sign-correlated pulses by 
10 taking into account parameters of spatial channel interaction. 

31. A method according to claim 30, wherein calculating the amplitudes of 
the sign-correlated pulses a single electrode model having exponential decays 
of the potentials at both sides of the electrode is used. 

15 

32. A method according to claim 31. wherein calculating amplitudes of the 
sign-correlated pulses, properties of a tri-diagonal matrix is used. 

33. A method of defining a channel specific sampling sequence having a 
20 pulse descriptive characterization, the channel specific sampling sequence 

being used to activate electrodes in a multichannel electrode array, each filter 
in a bank of filters associated with a channel having an electrode, the method 
comprising: 

i. deriving pulse amplitudes of the channel sampling sequence by 
25 sampling a signal waveform; 

ii. deriving the duration and number of pulses of the channel 
specific sampling sequence from a frequency associated with the charmers 
filter. 

30 34. A method according to claim 33, wherein the sampling is of a half 
period of a sinusoid between 0 and k. 
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35. A method according to claim 33, wherein the sampling is of a quarter 
period of a sinusoid between 0 and 7c/2 so that pulse amplitude distribution 
monotonically increases. 

5 36. A method according to claim 33, wherein the sampling uses biphasic 
current pulses. 

37. A method according to claim 33, wherein each filter is a bandpass 
filter. 

10 

38. A method according to claim 37, wherein the duration and number of 
pulses in the channel specific sampling sequence is derived from the center 
frequency of the charmers bandpass filter. 

15 39. A method according to claim 37, wherein the duration of the channel 
specific sampling sequence is one half of the period of the bandpass filter's 
center frequency. 

40. A method of deriving a weighting factor for a channel specific 
20 sampling sequence, the channel specific sampling sequence being used to 
activate electrodes in a multichannel electrode array, each filter in a bank of 
filters associated with a channel having an electrode, the method comprising: 
a. rectifying the output of each filter, creating a half-wave rectified 
signal; and 

25 b. determining a maximum amplitude of each half-wave in the 

half-wave rectified signal. 
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